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We study scalar perturbations to a Robertson- Walker cosmological metric in 
terms of a pseudo-Newtonian potential, which emerges naturally from the solution 
of the field equations. This potential is given in terms of a Green function for matter 
density fluctuations of arbitrary amplitude whose time and spatial dependence are 
assumed known. The results obtained span both the linearized and Newtonian limits, 
and do not explicitly depend on any kind of averaging procedure, but make the valid 
assumption that the global expansion rate is that of a Friedmann-Robertson- Walker 
model. In addition, we discuss the similarity to diffusive processes in the evolution 
of the potential, and possible applications. 
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I. INTRODUCTION 



Most observers would agree that the geometry of our universe is well- described by pertur- 
bations to a Robertson- Walker metric. These perturbations directly affect the information 
we receive from the universe outside of our own galaxy, and so a great deal of effort has 
been devoted to understanding their effects and evolution. 

Here we study the effects of a given matter distribution on the metric, and hence on the 
radiation that reaches us. Our goal is to find a way of expressing these perturbations directly 
in terms of the matter variables, in a way that may be useful for interpreting or modeling 
observations. Most simulations of observable effects — gravitational lensing, redshift, and so 
on — use some kind of radically simplifying assumption, such as nearest-neighbor Newtonian 
gravitation, weak density contrast (i.e. linearized theory), etc. Here we attempt to provide 
an expression which takes account of all cosmologically relevant effects, and which applies 
over a wide range of density contrasts. 

This paper expands on the results presented briefly in reference We show how to 
derive the Green function presented there for any value of the curvature parameter, as well 
as its reduction to the Newtonian limit and the similarity to simple diffusion. The principal 
results are summarized by equations ([11]), fl24]), (|5|), and (|8|) for the perturbed line element, 



pseudo-Newtonian potential, and Green function. 
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II. CONVENTIONS AND DEFINITIONS 



Units are chosen with G = c = 1 so that all quantities can be expressed as powers of a 
length. The metric is written as a perturbation on the static part of the Robertson- Walker 
background, 7 a b, with signature +2 and curvature parameter k G {—1,0, 1}: 



The coefficients and coordinates are dimensionless, so the scale factor a{rf) (which we write 
as a function of conformal time) is a length. Indices a,b,c, . . . run 0-3, except for i,j,...,n 
(the Fortran integers) which range over the spatial indices 1-3. An extremely useful source of 
formulas for the background metric and perturbations, in both conformal and proper time, 
is appendices A-D of Kodama and Sasaki, 1984 ||. Our notation is similar; in particular a 
prime (e.g. a') always means derivative with respect to conformal time. 

Naturally the perturbations h a b are assumed small compared to jab, which are of order 
unity. Of course this does not imply that the corresponding matter density fluctuations are 
small compared to the background. We use the scheme of Futamase 1988 H for parameteriz- 
ing the size of the perturbations and their derivatives, which is very much like the one used 
in the shortwave approximation for studying gravitational waves in a curved background 
(see Isaacson, 1968a ||J). With this scheme, orders of magnitude in our problem are: 



Formally, at least, perturbation of a curved background must parameterize the size of deriva- 
tives because the background radius of curvature is a natural length scale. In addition, many 
cosmological problems have a particle horizon length L. The parameter k represents the 
length scale of the perturbations relative to the particle horizon: k — l/L. In principle it 
can be either greater or less than one; in practice it is usually small. Specific restrictions 
will apply when we consider orders of magnitude more carefully in section |V|. 

Besides taking e 2 -C 1, it is also assumed that e 2 //t <C 1. As shown below, this implies 
that the matter inhomogeneities always move slowly (non-relativistically), and also that the 
effective stress-energy of the metric perturbations is small. 




(1) 




lab 



Oil) 



(2a) 
(2b) 
(2c) 
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III. THE PERTURBATION EQUATIONS 

The homogeneity and isotropy of the background 7^ allows separation of space and 
time dependence in the field equations, so it is useful to write the perturbations not as 
spacetime functions h a b but as a harmonic decomposition. (See, e.g., Harrison, 1967 [Q and 
Sasaki, 1987 ||.) The perturbations' spatial dependence is expanded in eigenf unctions, or 
normal modes, of the covariant Laplacian VjV* on the 3D static background 7^, reducing 
the gravitational field equations to a set of equations for the time- dependent amplitudes of 
the modes. The Appendix gives some general references to the (extensive) literature on the 
eigenfunctions. 

In this paper we focus on the scalar modes, representing perturbations to the metric and 
matter variables that can be written entirely in terms of solutions to the equations 



ViV'Qix , q) + q 2 Q = (3a) 

k = +l : q 2 = 0, 3, 8, . . . (n 2 - 1), (3b) 

k = : q 2 > 0, (3c) 

k = -1 : q 2 > 1 . (3d) 



This describes the sort of fields resulting from density fluctuations. Spatial derivatives are 
longitudinal with order of magnitude ViQ = 0(q) = 0(k~ 1 ). Vector (rotational) and tensor 
(transverse) modes are not included, since it is reasonable to suppose that most observables 
(such as lensing and redshift) are dominated by scalar effects. However, inflation-induced 
gravitational waves could have a significant effect on the anisotropy of the cosmic microwave 
background; see, e.g. Davis et al., 1992 Q. 

Representing the perturbations this way both simplifies and complicates the interpreta- 
tion of k. It simplifies, in that for a given scalar mode Q, k is just the reduced wavelength: 
k = q~ l . It complicates, in that for a real condensed object like a galaxy or cluster, k must 
be thought of as characterizing those wavelengths which contribute the largest physical ef- 
fect at a point. With this in mind, the expansion of the metric perturbations comes from 
Bardeen 1980 [|: 

27oo Jdfi(q) Q(x, q)A(rj, q), (4a) 
Jdfi(q) q^ViQix, q)B( V , q), (4b) 
27y Jdfi{q) Q{x,q)[H(r],q) + -H T (r],q)} 

+2 fdfi(q) q- 2 ViVjQ(x, q)H T (r), q) . (4c) 

Integrals over the measure n(q) stand for whatever operation expresses the completeness of 
the Q's, depending on whether the eigenvalue spectra in equations (|TJ) give a discrete or 
continuous set of modes. Most linearized calculations omit the integrals as being implicit, 
but in this paper we write them explicitly for clarity in what follows. 

Choosing longitudinal gauge, B = H T = 0, simplifies these expressions greatly. With 
only scalar perturbations this gauge fully specifies the metric and lets us think of the per- 
turbations as local length contractions and time dilations, making for easy comparison with 
Newtonian theory because it leaves the metric diagonal: 



^00 — 
hoi = 
hij = 



4 



h 



oo 



-2 JdfiQA, (5a) 
hoi = 0, (5b) 
ft^j = +2%j J dfx QH . (5c) 

Note that the amplitudes A(r], q) and H(rj, q) are supposed to have the same order of 
magnitude (e 2 ) as the metric perturbations themselves. 

(Other common choices for the gauge include harmonic and synchronous. Harmonic 
gauge is the usual choice for gravitational waves, and eliminates many higher-order per- 
turbation terms in the Einstein tensor, but the form of the metric is more complicated. 
Synchronous gauge is also common, and has the advantage of leaving the time coordinate 
unchanged. But it suffers from being under-specified, even when there are only scalar modes, 
allowing spurious gauge-mode solutions which can be difficult and annoying to identify and 
remove.) 

In longitudinal gauge the Einstein tensor is 



Coo 



3[(-) 2 + k] + 2 fdfi Q[q 2 H + 3-H' + 3k{A - H)], (6a) 
2jdfiViQ[^A~H'}, (6b) 



Ga = [(-) 2 - 2(-) - k] ltJ 



-2 7y fdy. Q \\q\A + H) + [(-) 2 - 2(-)](A - H) - —{A 1 - 2H>) + H"\ 

J I Zi (X (X Qj I 

- Jd^\7iS7 3 Q{A + H) . (6c) 

This includes all terms linear in h a b and its derivatives. It does not include nonlinear terms, 
which are C(e 4 ), 0(e 4 / K ), 0(e 4 /n 2 ), or smaller, the so-called pseudotensor terms discussed 
at the end of section 0. The stress-energy tensor for the matter is constructed by 
defining variables in the matter rest frame and then performing a Lorentz boost into the 
coordinate frame |J. The perfect-fluid background model is given (scalar) perturbations to 
density, pressure, and velocity as: 

P(77, x) = P(t]) + p{n) jdn(q) Q(x, q)U( V: q), (7b) 
Vi(ij, x ) = - ft dfi(q) q~ x ViQ{x, q)v(rj, q) . (7c) 

Actually, it is better to call A a density fluctuation because, as mentioned earlier, it is not 
necessarily small. However, the changes in velocity and pressure are always small, as shown 
below. (Note that the definition of the pressure perturbation n differs from Bardeen ||; 
this definition is easier to interpret when P = 0.) 

To first order in velocity, the boost to the coordinate frame gives components of T a b as: 
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Too = a p 
T oi = a 2 p 



1 + Jdp, Q(2A + A) 

(1 + er) Jdp q~ 1 V i Qv + Jdp, q' 1 W l Qv Jdp,QA 
a + JdfiQ(2cxH + U) 



(8a) 
(8b) 
(8c) 



(cr is the background pressure/density ratio.) Remainders are 0(v 2 ); the product term in 
the time-space component is kept because A may be greater than one. 

Most approaches to perturbation theory equate terms with equal orders of magnitude. 
But deferring order of magnitude arguments lets us exploit the harmonic decomposition of 
the field equations first, in essentially the same way as isolating coefficients in a Fourier 
series [see equation dl3|) 1. This has the advantage that the results obtained remain valid as 
the relative order of magnitude of terms in the field equations change — say as the density 
contrast is either diffuse (large k) or condensed (small k), or as its amplitude becomes large 
or small. 

Removing the trace from the Einstein equation Gij — 87rTjj = gives a result which looks 
familiar from linearized gravity (on flat spacetimes): 



2 J dpi {ViVjQ + ^q 2 HiQ){A + H) = 



which can only be true for arbitrary amplitudes A and H if 

A(r ] ,q) = -H(r ] ,q) + 0(e 4 ) . 



(9) 



(10) 



Using this in equations fllllfl and (|T|) for the perturbations and the line element puts the 
metric in the linearized pseudo-Newtonian form widely used for studies of gravitational 
lensing, etc: 



ds 2 
(f)(r],x) 



a 2 \-(l 



2<j))dr] 2 + (1 - 2cf)) lij dx i dx j } , 
l -h m = - Jdp(q) Q(x, q)H( V , q) + 0(e 4 



(11) 
(12) 



To obtain equations for H(rj, q), and thus <j)(rj, x), in terms of the matter variables we use 
equation ( ]T0| ) in the components of the orthogonality equation 



JdV Q*(G ab -8nT ab ) = 



(13) 



where dV is the proper volume element in the static 3-space 7^ (see the Appendix). The 
density fluctuations are governed primarily by the time component a = b = 0. Using 
equations (| 



and (lU) f° r the Einstein and stress-energy tensors gives: 



3(-)#' + {q 2 + 87ca 2 p - 6k)H = 4vra 2 pA . 

(X 



(14) 



To this level of approximation, the scale factor obeys the usual FRW equation for the 
background density p{rj) [compare equation (|20|) below]: 
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(15) 



In a formal sense, a must obey this equation in order that the integration in equation ( |13| ) 
gives a Dirac delta function in position when applied to the 0(1) (background) terms in the 
field equations, owing to the assumed completeness of the Q's. This allows us to make free 
use of substitutions from the background model for terms involving the scale factor and p 
in the following sections. 

Equations relating the metric perturbations to velocity and pressure perturbations follow 
similarly. The equations for II and v come most readily from the spatial and space-time 
components of ([13|) respectively: 

H" + 3(—)H' - (8na 2 pa + 2k)H = -4ira 2 pTl , (16) 
a 

q[H' + {-)H) = -4vra 2 p(l + <r)v[l + 0(A)] . (17) 

(Strictly speaking, the last equation holds only for q ^ — no real restriction, because q = 
is a constant mode, merely representing an improper definition of background.) We leave the 
A correction in ( |TT| ) inexplicit since products of eigenfunctions add nontrivial complications 
to the formalism, but this suffices for order of magnitude assessment. 
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IV. ORDERS OF MAGNITUDE 



Comparing the perturbations and their relative effects on the metric requires careful 
consideration of the size of the time derivatives. These are of order one in systems that 
are not gravitationally bound, H' = 0(e 2 ). However in bound systems we expect H' = 
0(e 3 1 ' k) — for instance, one can imagine an observer stationed a fixed distance R from the 
center of a collapsing dust cloud; changes in the potential are ~ (fi(V/R), where V is the 
speed of infall for the dust, or ~ e 2 (e/n) by order of magnitude. 



In equation fll4]) , the size of the density contrast depends on the ratio e/n [remember 
that q = The allowed regimes are labeled as linear or nonlinear density, depending 

on the size of A: 



O(H') 


0(A), K < 1 


0(A), K> 1 


LDR: c/k < 1 


e 2 


e 2 /n 2 


e 2 


NLDR: e/n » 1 


€ 3 /K 


e 2 /K 2 


not allowed (e> 1) 



Strong density fluctuations with large (super-horizon) scales are not allowed because they 
create potentials that cannot be treated as perturbations on the background metric. (So the 
results derived here, in particular the Green function given later, should not be expected to 
work in a strongly "tilted" universe.) However, the small-scale density fluctuations (NLDR 
with k < 1) should have an order of magnitude consistent with Newtonian theory. The 
prediction from the Poisson equation 



1 



2a 2 



oo 



(18) 



is also that A should be 0(e 2 /K 2 ). This holds for weak density fluctuations as well, again 
provided that the scale is small, since Newtonian physics can be expected to apply in a 
sufficiently small region of space. 

Orders of magnitude for pressure and velocity perturbations come from equations ( [TO] ) 
and (H3): 





0(H) 


0(v) 


LDR: e/n < 1 


e 2 


e 2 /n 


NLDR: e/n > 1 


e 4 /K 2 


€ 



These agree nicely with the simple Newtonian argument that for bound systems (NLDR) 
the velocity should be proportional to the square root of the potential, v ~ H 1 ' 2 ~ e, and 
the pressure should be n ~ Av 2 ~ e 4 //t 2 . More generally we can use the Euler equation for 
the hydrodynamics of a perfect fluid in a gravitational field: 



dv 



d t v + (v ■ V)v 



-p~ x VP 



V0 . 



(19) 



In the linear density regime the partial time derivative is most important, and comparing 
orders of magnitude we find v = 0(e 2 /k) and n = 0(e 2 ), as in the table above. In the 



S 



non- linear density regime p is of order e 2 /^ 2 , and again the estimates agree with the results 
from equations (0) and ([171) . 

The tables show, not surprisingly, that in any allowed (subhorizon) regime, the pressure 
and velocity perturbations are much weaker than the density fluctuations. So under these 
conditions the metric perturbations H(rj,q) are determined primarily by A(rj,q); that is, 
hydrodynamically the density fluctuations can be treated as the source. 

We also have to consider effects on the scale factor a, since it makes an implicit contribu- 
tion to any order of magnitude arguments. When nonlinear terms are kept in the Einstein 
tensor, their spatial average can be thought of as an energy density and used to construct an 
effective stress-energy tensor for the metric perturbations. (See, e.g., Isaacson, 1968b ||.) 
Dropping for a moment the requirement that a(r]) come from the background, on physical 
grounds we expect 

a(r l )=a FRW [l + 0(<e 4 /K 2 >)] (20) 

but clearly, the average over even a relatively small volume must be less than the maximum 
value; i.e. 0(< e 4 //t 2 >) <C c 4 /k 2 <C 1. So using the background scale factor does not alter 
the arguments about A, IT, and v, simply because the corrections are even smaller than any 
of the terms discussed previously. 
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V. GREEN FUNCTION 



Having determined that the scalar metric perturbations are determined primarily by 
the density fluctuations, we can relate them directly by solving equation ([14]) for the mode 
amplitudes H(r), q) and using the result in (0) to get 0(r/, x). Equation (|i~4"D is of the Riccati 
type; the standard solution given in handbooks like Gradshteyn and Ryzhik [HIT] is 



where 



H(rj, q) = H(r) , q)E(r} , r),q)- / du I(u, q)E(u, 7], q) (21) 



E(u, V,Q) = ^l exp[-(g 2 - 3k)C(u, 77)], (22a) 

/(M ' 9l = "l(^) AKg_) ' (22b) 
C(M) = !f*fi • (22c) 



3 

All of this, along with the definitions 



H(vo, q) = - JdV(y) Q*(y, q)<f>( Vo , y), (23a) 
A(u,q) = JdV(y)Q*(y,q)A(u,y), (23b) 



is inserted into equation (JT2[) . The integrals are re-ordered to contract the mode sums as 
much as possible, to give a kernel against which the initial conditions «/>(r/o, y) and source 
A(u,y) are integrated over space. After some manipulation this gives: 



^(t?, x) = JdV(y) G(r} , 77, x, y)(j)(vo, V) 

_^ r du f dV{ y) G (u, V , x, y)A(u, y) + 0(e 4 ) (24) 

where 

G(u, V , x, y) = ^±e 3kC ^) fd^q) Q(x, q)Q*(y, q)e^ c ^ ; (25) 

which is the central result of this paper: a Green function for metric perturbations due to 
scalar density fluctuations in a Robertson- Walker background. Equation (p4[) has several 
properties expected from an expression for a pseudo-Newtonian potential; for instance an 
overdensity A(u,y) > decreases </> in its neighborhood. Other simple cases can be checked 
with the help of the following formulas for special arguments of the Green function, which 
come from the completeness and orthogonality relations, respectively: 

G(u,u,x,y) = 5(x,y), (26) 

fdV G(u, V , x, y) = ^l e 3kC ^ . (27) 
J 1 a{ri) 
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The first shows that matches its initial conditions as rj — > i]q in fl24]) . The second, which 
holds as long as Q = const is an allowed mode, shows that is a function of time alone if A 
is, which means only that the background density has been shifted: p — > pA(i]); of course 
A = is also a solution. (Any spatially homogeneous form for the density fluctuation leads 
to a standard Robertson- Walker metric after transformation of the scale factor and time 
coordinate.) 

Specific forms for the Green function ( |2"5] ) come from choosing a particular representation 
(a coordinate system, in other words) for the Q's, given k. Regardless of the representation, 
in closed or open spaces (k = ±1) the integrals can be very ugly, especially for k = — 1. 
But for angles smaller than the curvature scale k = is a good approximation, and under 
inflation it would hold generally. In this case we can represent the Q's as plane waves and 
carry out the integrals explicitly. Using rectangular coordinates for q, the integral in ( |25| ) is 
separable, and completing the square in each exponent gives 

G k=0 (u,rj,x,y) = ^p^^xp 

The rest of the paper is devoted primarily to examining this formula and its use in equation 
(^4]) for the potential. Henceforth k is implicitly zero, unless noted otherwise. 



I — * -M 9 

\y — x\ 



(28) 
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VI. ANALOGY WITH DIFFUSION 



The most striking thing about equation ( 28) is that it looks very much like the Green 
function for diffusion in a uniform medium, and more generally, the displacement probability 



distribution for an isotropic random walk |TT|. For instance, we can compare the solution 
to the equation of heat conduction: 



dT 
~dt 



kV 2 T + c- 1 q . 



(29) 



(In this section, k and q represent thermal diffusivity and the rate of heat production per 
unit mass, not the reduced wavelength and mode numbers used in previous sections.) In 
an infinite uniform medium with no sources (q = 0) the solution for an arbitrary initial 



temperature distribution is [12 



T(f,t) 
G(f-f',t-t Q ) 



d 3 r' G{r 
1 



f',t-t )T(f',t ), 



[4tTk(* - t )] 3/2 



exp 



f\2 



An(t - t c 



(30a) 
(30b) 



Since there are no sources, the gravitational analog is the integral over initial conditions 
in equation (23) for the potential, using the Green function (^). If the integral for the 
potential is written with a physical volume element a 3 (r]o)dV, like the temperature integral 
above, we compare a~ 3 (r] Q )G(r] , rj, x, y) with equation ( |30b| ). Taking r] and rj close together 
to minimize the kinematic effects of expansion, the analog of the diffusivity and time factors 
is just 



K (t _ t(j ) ^ a 2 (??o)C(?7o, V) - (3^o) _1 At 



(31) 



where Hq is the Hubble parameter, evaluated at rjo in this particular instance, and At is 
the proper (rather than conformal) time interval corresponding to (77 — r] ) . This analog of 
the diffusivity is consistent with an alternate approach that comes from converting equation 
(|i"4]) for the mode amplitudes H(r],q) to an equation in (r],x), to form the analog of (p9|). 
Summed over modes and converted to physical variables (proper time, standard Laplacian, 
etc) it can be written as 



d t (a4>) = (ZH^V 2 



1 



Ho(aA) 



(32) 



where the Hubble parameter is taken to be a function of time. This is a diffusion equation 
for the potential, with time-dependent diffusivity and "specific heat capacity:" 



T 

K 

q 



—a<p , 

2H Q l , 
aA , 



(33a) 
(33b) 
(33c) 
(33d) 



in agreement with (|3l|). In simple kinetic theory, the diffusivity would imply a mean free 
path of Hq 1 for particles of average speed c = 1. In this approximation, changes in 
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travel at all speeds, as can be seen both from the analogy with diffusion and from the fact 
that equation (24) assigns a non-zero value to the potential everywhere, even for (spatially) 
localized density fluctuations. This superficial causality problem comes from dropping terms 
(9(e 4 //t 2 ) while constructing fll4l) . Ordinary diffusion cures the causality problem by making 
the diffusivity and specific heat capacity temperature-dependent |Tj|. In the gravitational 
case this means making Hq = a'a~ 2 a function of the metric perturbations — which it is, 
in an exact treatment, since the scale factor is affected by the spatial average of (V0) 2 
(cf equation (pOP). But causality does not represent a problem in the application of (p4|). 
however, because far from the source, "errors" in are extremely small — below the level of 
approximation for the calculation. 
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VII. NEWTONIAN LIMIT 



The diffusion analog is easiest to see when considering the role of the initial conditions. 
The Newtonian limit treats the opposite situation, where we consider compact sources (grav- 
itationally bound systems, for instance) evolving slowly and with negligible initial conditions. 
If the time derivative is ignored, the gravitational diffusion equation fl32|) becomes the Pois- 
son equation (18) for the potential, after replacing Hq by 8irGp/3. This makes it reasonable 
to suppose that the Green-function expression ( P4j ) for the potential, which is essentially de- 
rived from the gravitational diffusion equation, can be reduced to a Newtonian form under 
appropriate conditions, providing a useful check of the formula. 

We use the formulas for k = 0, but in fact the Newtonian limit holds in a sufficiently 
small region of any cosmological metric ||14|| , and so the results can be expected to hold for 
k = ±1 as well. (Even when k ^ the diffusion equation fl32|) still reduces to the Poisson 
equation when only the NLDR leading terms are kept.) Using the background equations for 
the scale factor and density we write the potential as: 

= - fdV \ \ Pdu —GA - G (f) 1 (34) 



'A i 



(I 



by analogy with the Newtonian form. The contribution of the initial conditions Go</>o is 
negligible provided that rjQ is taken far in the past and we assume that the early universe 
was smooth. Then to evaluate the time integral, we can imagine that the potential is to be 
measured at a point outside a nearby "lump" of matter. In this situation it is physically 
reasonable to suppose that the Green function will peak for values of conformal time u close 
to rj, corresponding to the time it takes for changes in the source to influence the observer. 
For small values of reversed conformal time w = i] — u the Green function is 

/ o A 3 / 2 

G= — - w~ 3/2 exp(-D /w) + O(w- 1/2 ), (35a) 
\ 47r a I 



3a' \y — x 



2 



This shows a peak at approximately w p = (2/3)D , with a width of about (2/5)D Q . Both 
are small provided that the matter distribution A(r],y) is such that the volume integral 
restricts D to small values. More precisely, an expansion in w is self-consistent if the width 
of the peak is small compared to the elapsed conformal time: D ^ (j] — i]q). With the flat 
background, and using r] = 0, this means that (omitting overall constants of order one) we 
must have 



\y- x 



2 

- < 1, or (36a) 



k 2 < 1 , (36b) 



where the last line follows from the "original" definition of n as perturbation size divided 
by particle horizon length. Thus the Newtonian limit holds for density perturbations that 
are highly localized, as one would expect. Furthermore, if the density field changes little in 
the time corresponding to the width of the peak, the source and expansion terms can be 
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evaluated at w = and moved outside the time integral. (Strictly speaking, w p might be a 
better choice, amounting to a notion of "retarded time," but the errors due to approximating 
the Green function make this academic.) Another change of variable to z = D /w puts the 
time integral in nearly standard form, and we find: 

/■rsn ! 

(37a) 




8=^- = lH /f «1. (37b) 
r] - 7] 4 a{r}- r] ) 

The z-integral is simply 7r 1 / 2 erfc(5 1//2 ), or, to a good approximation when 8 is small, just 
7T 1 / 2 . This along with the definition of Dq and the ever-present background equations for 
the scale factor reduce (|37a|) to the familiar Newtonian form: 




(38) 

To summarize, this formula holds when the initial conditions can be neglected, and the 
spatial dependence of the density contrast A(u, y) limits 8 to small values in equation (|37a|) , 



while the time dependence has a scale much larger than (any value of) D . In fact most grav- 
itationally bound systems satisfy these criteria quite well, and are also quite uninteresting 
from a cosmo logical standpoint. More interesting are situations where the time evolution of 
the density fluctuations makes a significant contribution to the metric. A theoretical descrip- 
tion requires extending the calculation just given to post-Newtonian order by expanding the 
Green function and source terms in a time-series. We hope to show in a forthcoming paper 
i"5fl how this is done, and under what conditions equation ( p4|) for the potential predicts 



significant deviations from the Newtonian (and LDR) approximations. 
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VIII. CONCLUSION 



In most situations of observational interest, the Green function expression for the 
pseudo- Newtonian potential offers a simple and relativistically correct way of calculating the 
metric perturbations, taking into account effects such as multiple (perhaps closely spaced) 
sources, deviations from the "thin lens" approximation, non-linear density evolution, and the 
cosmological expansion. The results can be applied in the calculation of observational effects 
such as lensing, redshift, and time-delay. We hope that this will stimulate exploration of 
situations that are difficult to treat with current techniques. A upcoming paper will examine 
the post-Newtonian limit of the Green-function expression, with attention to those situations 
which predict significant observable effects. 
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This appendix gives selected references to the literature on harmonic functions. Most 
calculations do not require an explicit representation; in spherical coordinates a schematic 
representation is: 



Q{x,q) = Iif\n,a)Y l 



lm\ 



x = (a, 6,<l)), 
q= (n,l,m), 

2 2 ; 

q = n — k, 

with the eigenvalue spectra in equations (|T|). The angular functions are ordinary spherical 
harmonics; good overviews of the properties of the radial functions and/or the Q's as a 
whole are in Harrison 1967 IBfl, Bardeen 1980 IBI, Kodama and Sasaki 1984 101, and Birrell 



and Davies 1989 [ IB] . The last writes the orthogonality and completeness relations as: 

fdVQ*(y,q)Q(y,p)=5(q,p) , 
Idfx Q*(x, q)Q{y, q) = S(x - y) , 



where dV = (t*) 1 ' 2 ^ is the proper volume element in the static background 3-space and 
d/u, is the measure associated with the eigenvalue spectrum. More detailed information, 
including explicit representations of the radial eigenf unctions and proofs of orthogonality 
and completeness can be found in (see also Harrison, above): Parker and Fulling 1974 ||17|| , 
and Abbott and Schaeffer 1986 II~ 



Finally, note that the literature aimed at problems in quantum field theory uses only 
scalar harmonics, while in general relativity an arbitrary tensor function may be composed 
of scalar, vector, and tensor harmonics. The early sections of Kodama and Sasaki || give a 
good explanation of the distinction. 



18 



REFERENCES 



* Electronic address: mark@relgyro.stanford.edu. 

t Electronic address: wagoner@leland.stanford.edu. 

* Electronic address: eric@bigbang.stanford.edu. 

[1] M. W. Jacobs, E. V. Linder, and R. V. Wagoner, Phys. Rev. D 45,R3292 (1992). 
[2] H. Kodama and M. Sasaki, Prog. T. Phys. Supp. 78,1 (1984). 

[3] T. Futamase, Phys. Rev. Lett. 61,2175 (1988); Mon. Not. R. Astron. Soc. 237,187 
(1989); T. Futamase and M. Sasaki, Phys. Rev. D 40,2502 (1989). 

[4] R. A. Isaacson, Phys. Rev. 166,1263 (1968). 

[5] E. R. Harrison, Rev. Mod. Phys. 39,862 (1967). 

[6] M. Sasaki, Mon. Not. R. Astron. Soc. 228,653 (1987). 

[7] R. L. Davis, et. al, Phys. Rev. Lett. 69,1859 (1992). 

[8] J. M. Bardeen, Phys. Rev. D 22,1882 (1980). 

[9] R. A. Isaacson, Phys. Rev. 166,1272 (1968). 
[10] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series, and Products (Academic 

Press, New York, 1983) pi 133. 
[11] S. Chandrasekhar, Rev. Mod. Phys. 15,1 (1943) Ch.l, §5. 

[12] A. L. Fetter and J. D. Walecka, Theoretical Mechanics of Particles and Continua 
(McGraw-Hill, New York, 1980) p415. 

[13] J. M. Hill, Solution of differential equations by means of one-parameter groups (Pit- 
man, Boston, 1982), Chapter 7; L. Dresner, Similarity solutions of nonlinear partial 
differential equations (Pitman, Boston, 1983), Chapter 4. 

[14] P. J. E. Peebles, The Large-Scale Structure of the Universe (Princeton Univ. Press, 
Princeton, New Jersey, 1980) pp37-41. 

[15] M. W. Jacobs and R. V. Wagoner (unpublished). 

[16] N. D. Birrell and P. C. W. Davies, Quantum Fields in Curved Space (Cambridge Univ. 

Press, Cambridge, England, 1989) pl21. 
[17] L. Parker and S. A. Fulling, Phys. Rev. D 9,341 (1974) appendix A. 
[18] L. F. Abbot and R. K. Schaeffer, Astrophys. J. 308,546 (1986) appendix. 



19 



